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ABSTRACT 

We analyze OGLE-2003-BLG-262, a relatively short, ?e = 12.5 ±0.1 day, microlensing event generated by 
a point-mass lens transiting the face of a K giant source in the Galactic bulge. We use the resulting finite- 
source effects to measure the angular Einstein radius, 6>e = 195 ± 17 //as, and so constrain the lens mass to 
the full-width half-maximum interval 0.08 < M/Mq < 0.54. The lens-source relative proper motion is fj, re i = 
27 ± 2kms _1 kpc -1 . Both values are typical of what is expected for lenses detected toward the bulge. Despite 
the short duration of the event, we detect marginal evidence for a "parallax asymmetry", but argue that this is 
more likely to be induced by acceleration of the source, a binary lens, or possibly by statistical fluctuations. 
Although OGLE-2003-BLG-262 is only the second published event to date in which the lens transits the source, 
such events will become more common with the new OGLE-III survey in place. We therefore give a detailed 
account of the analysis of this event to facilitate the study of future events of this type. 
Subject headings: gravitational lensing — stars: low-mass 



1. INTRODUCTION 

Immediately following the announcement of the first 
micro l ensing detectio ns (Alcocketal. 1993; Aubourg et al. 
Il993t lUdalski et alJ 119931) . three groups independently 
showed that one could measure the microlens angular Ein- 
stein radius, 



rel. 



AG 

c 2 AU 



mas 

'M 



77 ' (D 



from the deviations on the microlensing lightcurve in- 
duce d by the finite size of the source jGoulJ [T994at 
Nemiroff & Wickramasinghe 1994: IWitt & Maoll994l) . Here 
M is the mass of the lens and 7r re i is the lens-source rel- 
ative parallax. Although all three considered the case of 
a point-mass lens passing close to or over the face of the 
source star, the great majority of the actual Be measurements 
made over the ensuing decade used binary-l ens events in 
which the source passed over the binary c austic ( Albro wet alJ 
| 1999al l2000al I200U lAfonso et al.ll2000l lAlcock et alJl2000t 
lAn et al.ll2002l) . There has been only one single-lens event 
for which finite-source effects have yielded a measurement of 
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9e- This was the spectacular event MACHO-95-30, whose 
M4 III source of radius ~ 6Or was transited by the 
lens jAlcock et alJll997t) . In fact, of the more than 1000 
single-lens microlensing events discovered to date, only two 
have a measured 6>e by any technique. The other was the 
equally spectacular MACHO-LMC-5 whose source-lens rela- 
tive proper motion /i re i was measured by directly imaging and 
resolving t he source and the M-dwarf lens six years after the 
event ( Alcock et a fl l200l . The angular Einstein radius was 
then inferred from, 



% = Mrel fe, 



(2) 



where ?e is the Einstein crossing time, which had been mea- 
sured during the event. 

Measurements of 8e are important because they constrain 
the physical properties of the lens. For most events, the only 
measured parameter that is related to the physical properties 
of the lens is ?e, which (from eqs. Q and 0) is a combination 
of three such properties, M, 7r re i, and If Oe is measured, 
one then determines /i re i from equation (0, and the only re- 
maining ambiguity is between M and 7r le i (see eq. Q). In 
some cases, /j, e i is directly of interest. For example, measure- 
ment of the proper motion of the binary event MACHO-98- 
SMC-1 led to the conclusion that the lens was in the SMC 
itself rather than the Galactic halo dAfonso et al.12 000). 

In other cases, one can combine the measurement of Oe with 
other measurements or limits to further constrain the charac- 
ter of the lens. The most dramatic example of this would be 
measurement of the microlens parallax, 



7T E = 




(3) 



which can be determined eithe r by observing the e vent from a 
satellite in solar orbit (Refsdal 1966; Gould 1995) or from the 
distortion of the microlens lightcurve induced by the acceler- 
ated motion of the Earth jGoulJl9 92). If both 6e and 7r E are 
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measured, one completely solves the event. That is, 

Af=— 5-, Tr rel = 6 E n E . (4) 

K7T re l 

Unfortunately, while tt e has been measured for about a dozen 
events, only on e of these also has a firm measur ement of 9 E 
JAn et alJl20"02l) . although ISmith. Mao. Wozniakl d2003l> also 
obtained tentative measurements of both Be and ir E . 

Another type of constraint that can be combined with a 
measurement of 9 E is an upper limit on the lens flux, which 
can often be obtained from the lightcurve. This flux limit can 
be converted into a luminosity limit at each possible lens dis- 
tance. If the lens is assumed to be a main-sequence star, then 
using equation and some reasonable assumption about the 
source distance, one can p ut an upper limit on the lens mass 
(e.g., Albrow et al. 2000b). Even in the absence of any addi- 
tional constraints, equation Q can be combined with a Galac- 
tic model to make statistical statements about the lens proper- 
ties fe.g. JAlcock et alJl997t) . 

The principal reason that most 6 E measurements come 
from binary lenses and that single-lens measurements are ex- 
tremely rare is that the ratio p, of the angular source radius, 
#*, to the Einstein radius, 



is usually extremely small. At the distance of the Galac- 
tic bulge, even a clump giant has an angular radius 6>* ~ 
6/Ltas, and main-sequence stars are an order of magni- 
tude smaller. By contrast, typical Einstein radii are 8 E ~ 
31O//as[(M/O.3M )(7r re i/4O/ias)] 1 / 2 . Hence, the probability 
that the lens will pass directly over the source, which is what is 
required for substantial finite source effects (iGould & Welch! 
1996), is very small. By contrast, binary lenses, with their ex- 
tended caustic structures, have a much higher probability of 
generating finite-source effects. 

However, new microlensing surveys are beginning 
to alter this situation. In particular, the new phase of 
the Optical Gravitat ional Lens Experiment, OGLE-III 
( Udal ski etail l2002al) . with its dedicated 1.3 m tele- 
scope and new 35' x 35' field, 0."26 pixel, mosaic CCD 
camera and generally excellent image quality is gener- 
ating microlensing alerts at the rate of 500/season (as 
reported by the OGLE-III Early Warning System, EWS, 
http://www.astrouw.edu.pl/~ogle/ogle3/ews/ews.htmll, 
roughly an order of magnitude higher than previous sur- 
veys. This high event rate is itself enough to overcome 
the low, O(p), probability of a source-crossing event, and 
so to generate a few finite-source affected EWS alerts per 
year. Moreove r, because EWS relies on image-subtraction 
(Wozniak 2000), it is sensitive to extremely high magnifica- 
tion events of relatively faint sources, which have a higher 
chance of a source crossing than do typical events. 

OGLE-III is able to generate this high event rate only by 
reducing its visits to individual fields below 1/night. Hence, 
it would not customarily observe the alerted event during the 
lens transit of the source. However, several groups, includ- 
ing the Probing L ensing Anomalies NETwork (PLANET, 
Albro w et alJll998l) . the Microlensing Planet Search (MPS, 
Pvhie et al. 19991 and the Microlensing Follow-Up Network 
(/iFUN, http://www.astronomy.ohio-state.edu/~microfun/! 
intensively monitor alerts from EWS and also from the 
Microle nsing Observatio ns in Astrophysics collaboration 
(MOA, iBond et all 120011) . primarily to search for planets. 



High-magnificat ion events are t h e most s ensitive t o planetary 
pertur bations iGould & Loebl 119921 iGriest & Safizadehl 
1998), so these groups tend to focus on these events, partic- 
ularly their peaks. As a consequence, there is a good chance 
they will detect finite-source effects when they occur. More- 
over, OGLE-III diverts time from its regular field rotation 
(survey mode) to especially interesting events (follow-up 
mode) and so can itself also directly detect these effects. 

Here we report observations of EWS alert OGLE-2003- 
BLG-262, which exhibited clear indications of finite-source 
effects near its peak on 2003 July 19. By fitting this event to a 
single-lens finite-source model, we measure the 9 E and so p, K \. 
We use this information, combined with a measurement of t E 
to constrain the mass of the lens. We also present marginal 
(> 3 a) evidence for an asymmetry which, if due to parallax 
effects, would imply that the lens was a brown dwarf. How- 
ever, we argue that the observed asymmetry is either due to 
statistical fluctuations, a weak binary lens, or acceleration of 
the source. Our analysis provides a framework in which to an- 
alyze future finite-source single-lens event, which should be 
considerably more common due to the higher rate of alerted 
events. 

2. OBSERVATIONAL DATA 

The microlensing event OGLE-2003-BL G-262 was iden- 
tified by the OGLE-III EWS dUdalski et alJH99l on 2003 
June 26, i.e., more than three weeks before peak, which oc- 
curred on HJD' = HJD-2450000 = 2839.84 over the Pacific 
Ocean. OGLE-III observations were carried out with the 1 .3- 
m Warsaw telescope at the Las Campanas Observatory, Chile, 
which is operated by the Carnegie Institution of Washington. 
These comprise a total of 170 observations in / band, includ- 
ing 68 in the 2001 and 2002 seasons. The exposures were 
generally the standard 120 s, except for three special 40 s ex- 
posures on the peak and following night when the star was 
too bright for the standard exposure time. Photometry was 
obtained with the OGL E-III im age subtraction tech nique data 
pipelin e (lUdalski et al.ll20 02a) based in part on the Wozniak 
(2000) DIA implementation. The source had also been mon- 
itored by OGLE-II and was found to be very stable over four 
previous seasons (April 1997 - October 2000). 

Following the alert, the event was monitored by /iFUN from 
sites in Chile and Israel. The Chile observations were car- 
ried out at the 1.3m (ex-2MASS) telescope at Cerro Tololo 
InterAmerican Observatory, using the ANDICAM, which 
simultaneousl y ima ges at optical and infrared wavelengths 
(DePo v et alJ 12003). During the seven nights from HJD' 
2838.5 to 2844.8, there were a total of 45 observations in /, 
4 in V, and 28 in H. The / and V observations were gen- 
erally 300 s, although the exposures were shortened to 120 
s during the three nights from 2839.6 to 2841.8. The in- 
dividual H observations were 60 s and were grouped in 5 
dithered exposures, which were taken simultaneously with 
one 300 s V or / exposure or with two 120 s / exposures. 
All images were flat fielded using sky flats for V and /, and 
dome flats for H. P hotometry was obtained with DoPHOT 
JSchechter. Mateo & Sahal 19931) for all V, I, and H images. 

After reductions, the contiguous groups of 5 (or 10 in the 
case of back-to-back V and / exposures) H points were av- 
eraged into single data points to yield the above-stated 28 
points. 

The /iFUN Israel observations were carried out on the Wise 
lm telescope at Mitzpe Ramon, 200 km south of Tel- Aviv, 
roughly 105° east of Chile. During the nights of ~ 2839.3, ~ 
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FIG. 1. — Photometry of microlensing event OGLE-2003-BLG-262 near 
its peak on 2003 July 19.34 (HJD 2452839.83). Data points are in / (OGLE: 
empty circles; /iFUN Chile: empty triangles; /iFUN Israel: crosses), V 
(/iFUN Chile: filled squares; /iFUN Israel: filled circles), and H (/iFUN 
Chile: filled triangles). All bands are linearly rescaled so that F s and F/, are 
the same as the OGLE observations, which define the magnitude scale. When 
the lens is close to or inside (vertical lines) the source, the lightcurves are ex- 
pected to differ due to limb darkening (LD). The solid curve shows the best fit 
model for the /-band curve. The fact that the //-band points near the peak are 
below this curve is in qualitative accord with the lower LD in H compared to 
/. The dashed curve shows the lightcurve expected for the same lens model, 
but a point source. 



2841 .3, and ~ 2842.3, there were a total of four observations 
in / and three in V. The exposures (all 240 s) were obtained 
using the Wise Tektronix IK CCD camera. Data were flat- 
fielded and zero corrected in the usual way, and photometry 
obtained with DoPHOT. 

The position of the source is R.A. = 17 h 57 m 08?51, decl. = 
-30°20'05f'l (J2000) (l,b = 0.41918,-3.46935), and so was 
accessible for most of the nights near peak from Chile, but 
only a few hours from Israel. Unfortunately, due to a commu- 
nications mixup, /iFUN Chile observations on the peak night 
were bunched in a narrow time interval. Happily, when these 
are combined with the two OGLE observations and the one 
/-band and one V-band /iFUN Israel observations, the rising 
half of the peak is still clearly traced out. See Figure^ 

We emphasize that while three of the datasets have rela- 
tively few points, two of these small datasets actually play 
crucial roles. The three post-peak /iFUN Israel / points serve 
to align this dataset with the two larger / datasets and so en- 
able the first point (on 2839.38) to directly test the near-peak 
finite-source profile, which otherwise would be determined by 



a single compact set of points. See Figure^ The four /iFUN 
Chile V points allow determination of the color of the source 
and so permit one to estimate the source size and thus the 
proper motion and angular Einstein radius. See §|4] With only 
three points, two of which are nearly coincident, the /iFUN 
Israel V data do not contribute significantly to the fit because 
they are absorbed by two fitting parameters, F s and F\,. How- 
ever, they are included here for completeness. 

The source lies in one o f the OGLE-II cal- 
ibrated photometry fields (Udalski et al. 2002b) 
( ftp://bulge.princeton.edu/ogle/ogle2/maps/bulge/) and 
this allows us to place it on a calibrated color magnitude 
diagram. See Figure |2] The source lies on the red giant 
branch, about 1 mag brighter than the clump and about 0.2 
mag redder. It therefore has considerably larger angular 
radius than typical microlens sources, and this, together with 
the high magnification of the event, considerably increased 
t he chance fo r signifi cant finite source effects. 

Sumi et al. (2003) measured the proper motion of the 
source (relative to the frame of the Galactic bulge) and found 
(McwM*) = (°- 45 ± 0.41,-5.75 ± 0.40) mas yr" 1 . When 
corrected to the Tycho-2 frame, this becomes (/x a , s ,/x<5,j) = 
(-2.9,-12.4) masyr -1 . 

3. FORMALISM 

3.1. Finite-Source Effects 

In most cases, the lensed star is regarded as a point source 
because the angular size of the source is negligibly small com- 
pared to the angular separation of the source and the lens. The 
magnification is then given by ( Pac zvriskil 19861) . 



A(u) = 
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m(m 2 + 4) 1 /2' 
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FIG. 2. — Calibrated color-magnitude diagram (CMD) of a 10' square 
around OGLE-2003-BLG-262 taken from OGLE-II photometry well before 
the event. The source (marked with a black triangle) is about 1 mag brighter 
and 0.2 mag redder than the centroid of the clump giants (marked with a red 
circle). The fit shows negligible blending, so the apparent source position on 
the CMD is virtually identical to its true (deblended) position. 
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FIG. 3. — Finite source functions Bq(z), #1/2(2) andfii (z) given by eqs. 1101 . 
1191 and 1161 . For p <C 1, the limb-darkened magnification is very well repre- 
sented by Ai,j(n|p) = A(u)[Bq(z) — TB[(z)], where p is the source size and u is 
the lens-source separation, both in units of Be, T is the linear limb-darkening 
coefficient, and z = u/p. 



where u is the projected source-lens separation in units of 
the angular Einstein radius 9%. However, this approximation 
breaks down for w < p. Finite-source effects then dominate. 

If the source were of uniform brightness, the total mag- 
nification would simply be the mean magnification over the 
source, 

A uni (u\p) = Wo[(u/p)\p;A(x)], (7) 

where 

W„ \z\p;f(x)] = - [dd [ ir r( l-r 2 )"^ 2 f(py/r 2 + z 2 - 2rzcos 6] 
^ Jo Jo 

(8) 

iWitt & Mad (119941) gave an exact evalua tion of this exp res- 
sion in closed, albeit cumbersome, form. Gould (1994i ad- 
vocated a simple approximation to equation (|8j, 

A uni ( M |p) ~A{u)B {u/p), B (z) = zpWo , (9) 

which follows from the fact that A(u) ~ u~ l when u <C 1. Note 
that B o depends on p only through the ratio z = u/p. How- 
ever, |Gou!d|lll994a) did not explicitly evaluate Bo nor did he 
demonstrate the range of validity of the approximation i|9}. It 
is straightforward to show that, 



B {) {z)=-zE(k,z), 

7T 



i(z _1 ,l), 



(10) 



where E is the incomplete elliptic integr al of the second kind, 
and w here we follow the notation of Gradshtevn & Rvzhik 
cum. Using the expansion A(u) = u 1 [1 + (3/8)m 2 +. . .], and 
after some algebra, one may show that to second order in p, 



A uni (u\p) = A(u)B (z) 



z= -, 



(11) 



Q(z) = 



7-8z 2 -4(l 



E(k,z) 



(12) 



and where F is the incomplete elliptic integral of the first kind. 
We find numerically that -0.38 < Q(z) < 1, where the limits 
are saturated at z = 0.97 and z = respectively. For OGLE- 
2003-BLG-262, p 2 /8 < 5 x 10" 4 , which is about an order of 
magnitude smaller than our photometric errors. The zeroth- 
order approximation is therefore appropriate here and, we 
believe, is likely to be appropriate in most other cases as well. 

3.2. Limb Darkening 

Real stars are not uniform, but rather are limb-darkened. 
For simplicity and also because the quality of the data does 
not warrant a more sophisticated treatment, we adopt a one- 
parameter linear limb-darkening law for the surface bright- 
ness of the source, 



Sx(#) = Sx 



l-Tx 1 



-cos ■& 



(13) 



where $ is the angle between the normal to the stellar surface 
and the line of sight, S\ is the mean surface brightness of the 
source, and the T\ is the limb-darkening (LD) coefficient for 
a given wavelength band A. The factor 3/2 originates from our 
requirement that the total flux be F tot .\ = tt6 2 S\. 
The LD magnification is then (exactly), 

A lld (u\p) = Wo[(u/p)\p;A(x)] (14) 

-r {Wo [( U /p)\p;A(x)] - 1.5 Wi [( U /p)\ P ;A(x)] } . 

However, we adopt the same simplifying approximation as 
above and write, 



where 



A Ud ( M | i o)~A(M)[fio(z)-rBi(z)] ; 



B 1 (z) = B (z)--zpW 1 [z\p;x~ l ] 



(15) 



(16) 



Figure 13711 shows Bo, B\n (see below) and B\ as functions 
of z. Note that Bq(z) — > 1 and Bi(z) — » in the limit z — > oo so 
that the magnification (eq. 1151 ) reduces to the point-source 
case. In the opposite limit, z — > 0, equations dlOt and dl6l 
reduce to B (z) —> 2z and B\(z) -> (2-37r/4)z, so that A fin (0) = 
2/p [l +(37r/8— l)rl. Hence, the peak magnification depends 
primarily on p and only weakly on T. 

In high-precision LD measurements, it is generally ac- 
cepted that a two-parameter square-root LD la w is more ap- 
propriate to describe b rightness profiles of stars ( Albro w et alJ 
1999a; Fields et al. 2003) than equation Jl 31 although it is not 
used in the present work. Therefore, for completeness we ex- 
tend the above formalism to a two-parameter square-root LD 
law in the form of 



Sx(#) = Sx 



! - r v ( 1 - ^cos ■& j - A A (l - ^cos I/2 ■& 



(17) 

where A\ is the additional LD coefficient for a given wave- 
length band A. The magnification can then be approximated 

by, 

A sqrt i d (K|A> ~A(a) [B (z)-TB 1 (z)-AB l/2 (z)] , (18) 
where 

By 2 (z) =B (z)- 7 zpW 1/2 [z\p;x~ l ] . (19) 
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Table 1. OGLE-2003-BLG-262 Fit Parameters 



Parameter 


Free Fit 


Fixed LD 


Fixed LD 


& 7T E 


Value 


Error 


Value 


Error 


Value 


Error 


/nfdavsl 


2839.8411 


0.0015 


2839.8415 


0.0015 


2839.8424 


0.0014 


^0 


0.0365 


0.0005 


0.0362 


0.0004 


0.0360 


0.0004 


'b\ Utt ; 9 / 


12.5309 


0.0945 


12.5568 


0.0941 


12.6181 


0.0916 





0.0605 


0.0010 


0.0599 


0.0005 


0.0595 


0.0005 


r v 


0.8515 


0.2069 


0.7200 




0.7200 




Ti 


0.6118 


0.1499 


0.4400 




0.4400 




r« 


0.0975 


0.2028 


0.2600 




0.2600 




te.ii 


-0.8572 


0.3130 


-0.8335 


0.3120 


0.0000 




(F b /F s ) h 


-0.001 1 


0.0095 


0.0028 


0.0093 


0.0083 


0.0091 


(F b /F s ) h 


-0.0275 


0.0175 


-0.0134 


0.0172 


-0.0027 


0.0168 


(F b /F s ) h 


0.1865 


0.0783 


0.2283 


0.0746 


0.2361 


0.0749 


(F b /F s )v 2 


0.0122 


0.0481 


0.0192 


0.0478 


0.0296 


0.0479 


(F b /F s ) v ; 


0.0406 


0.1688 


0.1251 


0.1465 


0.1191 


0.1473 


(Fb/F^H 


-0.0048 


0.0176 


-0.0114 


0.0175 


-0.0009 


0.0172 


x 2 


233.50 




252.66 




259.76 





NOTE. — Observatories: l=OGLE, 2=^FUN Chile, 3=^FUN Israel 



3.3. Parallax Effects 

Microlensing events are fit to 

F{t) = F s A[u(t)]+F h , (20) 

where F s is the source flux, Ft is the blended background light, 
and 

"(0= Vlr(t)l 2 +W(t)] 2 . (21) 
Conventionally, rectilinear motion is assumed, 

r(f)=— , 0(0 = »<>. (22) 

Hence, the simplest fit has five parameters, F s , F h , the impact 
parameter uo, the time of closest approach to, and the Einstein 
timescale t E . However, even if the source and lens are in rec- 
tilinear motion, the Earth is not. Thus, strictly speaking one 
should write 

r(f)= t —^- + ir E uan(t),+ir E ,j_a ± (t), (23) 

/3(t) = u —K nill a 1 _(t) + Tr B ,xa {l (t). (24) 

Here a = (an,a±) is the difference in the Earth's position (pro- 
jected onto the plane of the sky and measured in AU) relative 
to what it would have been had the Earth maintained the ve- 
locity it had had at to, and the an direction is defined by the 
direction of the Earth's (projected) acceleration at to- 

Choosing the Earth frame at the peak of the event as the in- 
ertial frame is certainly not standard procedure. It is more 
common, and mathematically more convenient, to use the 
Sun's frame. However, for relatively short events t E < yr/27r, 
the parallax effect is quite weak, and it is only possible to 
measure one component of tt e = (7Te,|| , te,_l), namely the par- 
allax asymmetry, which is the component {ir E ii) of the paral- 
lax parallel to the Earth's projected acceleration at t he peak 
of the event (iGould. M iralda- Escude & Bahcalll 19941) . In this 
case, uo, to, and t E as seen from the Earth at the event peak 
are very well defined by the fit to the event without parallax, 
whereas these quantities as seen from the Sun are impossible 
to determine. For these short events, a± ~ 0, and the impact 
of a || through /3 is undetectable because it is absorbed into uq, 
t E , F s , and F b . Equations d23i and J24i then reduce to, 

T(0=^ + 7r E ,||a||(0, Pit) = "o- (25) 



4. MODEL FITTING 

We begin by fitting the event taking account of both LD and 
parallax. There are then a total of 20 free parameters: 12 pa- 
rameters for F s and F/, from each of the six observatory/filter 
combinations, 3 LD parameters, one each for /, V and H, the 
basic microlensing parameters to, uq, and t E , as well as the 
source size, p, and the parallel component of the parallax, 
7r E || . We consider the possibility of a correction for seeing, 
but find no correlation of the residuals of the /iFUN Chile / or 
H data with seeing. The source is quite bright (see Fig.|2j and 
it is virtually unblended (see below), so it is quite plausible 
that there would be no seeing correlations. We set a minimum 
error of 0.003 magnitudes for all observations, regardless of 
what value the photometry programs report. We then rescale 
the errors for the OGLE, /iFUN Chile /, and H by factors 
of 1.62, 1.12, and 1.83 respectively, in order to force x 2 /dof 
to unity. There are too few points in each of the remaining 
observatory/filter combinations to permit accurate rescalings, 
and the actual total x 2 values for these are consistent with the 
reported errors being correct. 

We minimize \ 2 using Newton's method (e.g. JPress et al.l 
which guarantees that one has found a local min- 
imum because the derivative of \ 2 with respect to each 
parameter is zero. In contras t to caustic-cr ossing bi- 
nary lenses ( Albrow et al.| [1999bt iDo minikl I1998I) . and to 
space-based JGouldll l994b: Refsdal 1966) and ground-based 
( Smith, Mao & Paczvriski 120031) parallax measurements for 
which there can be multiple local minima, standard mi- 
crolensing (even when modified by inclusion of finite-source 
effects) is expected to have a single global minimum. We 
nevertheless checked for multiple minima by adopting sev- 
eral initial trial solutions that were consistent with point- 
source/point-lens fits to a data set that excluded the peak. All 
converged to the same solution. 

We initially allow the three LD coefficients to be free pa- 
rameters. We find fit values and errors (Tv,Tj,Th) = (0.85 ± 
0.21 , 0.6 1 ± 0. 15 , 0. 10 ± 0.20) (see Tab. Q). These errors are 
all relatively large. The values therefore appear only mildly 
inconsistent with those of EROS-BL G-2000-5, (IV, T/, r H ) = 
(0.72.0.44,0.26) jFields et alJl2003l) . a slightly redder source 
with much better measured LD. It is then somewhat shocking 
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FIG. 4. — Model-independent color changes due to limb-darkening. A 
linear regression of H on / flux is performed at high z (z > 1.7) to put the two 
passbands on the same scale and to remove the small blending difference. 
Then /—//is measured at each point and the measurements for each day are 
averaged, except for HJD' ~ 2840.5, (z ~ 1.25), which is broken into two 
bins. The curve is 0.5Z?i(z), which is the expected form of this magnitude 
difference for a linear limb-darkening difference Vj—Vh = 0.5. 



rescaled to have an F s and Ft, equal to those of the OGLE 
data set. The three / band data streams should then follow the 
same lightcurve, whose best fit model is shown by the solid 
curve. However, the V and H band points should deviate from 
this curve during the source crossing, \t—to\ <pt e = 0.76 days, 
because of LD. There are not enough data in the V band to test 
this. As mentioned above, the H band cluster of points near 
the peak clearly lies below the curve, probably by too much. 
We measure a parallax asymmetry, 

7T E .|| =0.83 ±0.31. (27) 

That is, parallax is formally detected at the 3 a level. More 
specifically, A\ 2 = 7 relative to enforcing tt e » = 0. To illus- 
trate the strength (or lack thereof) of this detection we show 
in Figure |5]the fit to the data enforcing tt e ii = 0. The lower 
panel of this figure shows the residuals together with their ex- 
pected form for ir E u = 1 . When we first constructed this figure 
on about HJD' = 2870, we realized that there might still be 
time to test the reality of this parallax detection. OGLE ob- 
servations were then intensified from one every several days 
to one or two per day. These additional observations did not 
tend to confirm the detection, but also did not firmly con- 
tradict it. Hence, the parallax detection remains ambiguous. 
All previous events with firm parallax detections had Einstein 
timescales at least 5 times longer than this one, so it would 
have been remarkable if we had obtained a robust detection. 
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to discover that there is a net penalty of A\ 2 =19 for enforc- 
ing the EROS-BLG-2000-5 LD values. A major part of the 
problem is that while the errors in the individual LD parame- 
ters are large, the data strongly demand a large LD difference 
AT = T,- T H = 0.5 1 ± 0.09 when T 7 is held fixed at 0.44. That 
is, although the errors on the individual determinations of Tj 
and Th are large, they are strongly correlated, such that the 
difference Ar is much better determined. This in turn can be 
traced to the fact that there is a color offset A{I-H) = -0.03 
at the peak, which is clearly visible in Figure^and which the 
fitting routine ascribes to the source having much more LD in 
/ than H and hence being relatively blue in the center. See 
Figure |4] However, since the measurement seems to contra- 
dict what is otherwise known about LD, and derives primarily 
from a single cluster of data points, which may be subject to 
common systematic error, we choose to fix the three T's at the 
above stated EROS-BLG-2000-5 values. We thereby lose any 
independent LD information. This is not a major loss since 
our errors are too large to be competitive with other measure- 
ments (e.g. Fields et al. 2003). Our main concern is that what- 
ever problem may be corrupting the LD could also impact the 
measurement of the parameters that we are most interested in 
measuring, which are principally p and t E . In fact, by enforc- 
ing these T's, p is changed by only 1.6% and t E by only 0.3%. 
Since enforcing the LD parameters has no practical conse- 
quences (other than the loss of LD information), we adopt the 
EROS-BLG-2000-5 value. 
We then find, 

p = 0.0599 ±0.0005, f E = 12.557 ±0.094days. (26) 

Figure [2 shows the fit to the data in the region of the peak. 
All six observatory/band combinations have been linearly 




FIG. 5 . — Similar to Fig.0 but now a full view of the lightcurve of OGLE- 
2003-BLG-262 over about 5 Einstein timescales. The fit does not include 
parallax, and the residuals (lower panel) show an asymmetry such as would 
be induced by acceleration of the Earth parallel to the direction of lens motion 
(solid curve). 
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Moreover, as we discuss in §[6] the sign of the effect is oppo- 
site to what would be produced by the expected lens-source 
kinematics, while the effect itself could be produced by xal- 
larap or by lens binarity. 

The errors shown in Tabled are \fcu where c, ; is the /y'-th 
element of covariant matrix, and the correlation coefficients 



defined as c,- , 
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md(F b ) h . As 



where parameters are fo, uq, fg, P, te,||> (Fs) 
expected from experience with standard microlensing events, 
F s and Ft are extremely correlated, and these are both highly 
correlated with uq and fe. What is new in equation d28i is 
that p is also highly correlated with these other four parame- 
ters. The fundamental reason for this is that all five of these 
parameters are symmetric in (t— to). By contrast ir E n is only 
weakly correlated with the other parameters. 

5. CONSTRAINTS ON THE EVENT 
5.1. Angular Einstein Radius 0e 

As discussed by Albrow et al. (2000a), one can determine 
6* from the source's dereddened color and magnitude [(V — 
I)oJo\s by first transforming fr om (V - 7)n to ( V - K) p us- 
ing the color-color relations of Bessell & Brett ( 1988), and 
then applying the empirical relati on betwee n color and sur- 
face brightness to ob tain 9* (Ivan Belle! 19 99). 

Again following Albrow et al. (2000a), we determine [(V — 
I)o,Io]s from the measured offset of the unamplified source 
(as determined from the microlensing fit) relative to the cen- 
troid of the clump giants on an instrumental CMD, the latter's 
dereddened color and magnitude being regarded as "known". 
We measure this offset to be 

A/ = /-/ c i ump = -1.06, A(V-/) = (V-/).-(V-/) c i urap = 0.15. 

(29) 

In general, the source may be blended, so that the V and / of 
the source derived from the microlensing fit will not necessar- 
ily agree with those of the object identified as the "source" in 
an image taken at baseline. Hence, one cannot in general de- 
rive the offset from a CMD constructed from such a baseline 
image. However, in this case, there is essentially no blending, 
so the offset in the baseline calibrated CMD shown in Fig- 
ure |2]is virtually identical (within 0.02 mag) to that given in 
equation (I29> . 

The "known" values of [(V-/)o,/o]ciump have recently come 
under dispute. The basic problem is that the previous calibra- 
tions of these quantities relied on a number of steps, in each of 
which it was assumed that the ratio of total to selective extinc- 
tion was Ryj = Ay /E(V - 1) ~ 2.5. However, using this same 
value. iPaczvnskil d 1 9961 ) and lStutz. Pooowski & Gouldl(ll999l) 
found respectively that the colors of bulge clump giants and 
RR Lyrae stars w ere anomalous relative to local populations. 
Popowski (2000) then proposed that these anomalies could be 
resolved if the dust toward this line of sight were itself anoma- 
lous, with Ryi ~ 2.1. Udalski (2003) then demonstrated that 
this was very likely the case based on OGLE-II data. While it 
would be both worthwhile and feasible to retrace all the steps 




-2-10 1 
Log (M/M ) 

FIG. 6.— Constraints on the lens mass of OGLE-2003-BLG-262. The 
curves show the relative probability of different lens masses given the mea- 
surement of $e = 195 /ms and the mass distribution along the line of sight 
as predicted by the Han & Gould 1 1 995 , 2003 ) model. The solid and dashed 
curves show the probability for bulge-bulge and disk-bulge combinations of 
lenses and sources. The bold curve is their sum. The constraints arising from 
the determination of /i re i (equivalently fg) would be extremely weak and are 
not incorporated here. 



that led to the old calibration in light of this revised Ryi, the 
magnitude of this project lies well beyond the scope of the 
present work. Pending such a revision, we adopt a simpler 
approach. 

The distance to the Gal actic cen t er has now been measured 
geometrically bv lEisenhauer et alJ 12003) to b e Rq = 8.0 ±0.4 
kpc ba sed on the "visual-binary" method of Salim & Gould 
( 1999). Bulge stars are of similar metallicity to local stars, so 
the clump should be of similar color to the Hipparcos clump 
stars (V-I)o ~ 1 .00. (Recall that it was the apparent failure of 
this expectation that led to the discovery of anomalous extinc- 
tion.) The /-band luminosity of clump stars does not depend 
strongly on age (until the stars are so young that the turnoff 
luminosity approaches that of the horizontal branch). Hence, 
the bulge clump stars should have approximately the same 
Mi as the Hipparc os sample. For this we ad opt M/ = -0.20, 
the value found by Paczvhski & Stanekl ( 119961) for their 70 pc 
sample (and prior to their reddening correction which we con- 
sider to be substantially too large.) Hence, in lieu of a more 
thoroughgoing calibration, we adopt 



t(V-/)o,/o]cium P = (LOO, 14.32). 



(30) 



Combinin g equ ations \29\ and (I30i and applying the 
Belld 3l999) relation, we find 



van 



= 11.7 ± 1.0 fjas, 



(31) 



where th e error comes pri marily from the 8.7% intrinsic scat- 
ter in the van Belle ( 1999) relation. 

This evaluation would appear to depend on the assumption 
that the source suffers exactly as much extinction as a typical 
clump star, which it might not, either due to highly variable 
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extinction or to the source lying well in the foreground and 
so in front of a large fraction of the dust. In fact, if it were 
determined that the extinction toward the source were greater 
than to the clump by AE(V-I) = 0.2 or less by AE(V-I) = 
-0.6, the estimate of 6>* would change less than 3%. This 
is because the changes in the inferred surface brightness and 
luminosity lead to changes in the source-size estimate that go 
in opposite directions. 
Combining equations \26\ and J3 II . we obtain, 

9 E = 195±17^as, (32) 
p lel = 5.63 ± 0.49 masyr" 1 = 26.7 ± 2.3kms _1 kpc" 1 . 

We now use this measurement of 8e, in conjunction with its 
definition, equation 0, to write the source-lens relative par- 
allax as a function of the lens mass, 



7r re i(A0 = ^=4.8Atasl ^~ 
kM V M Q 



(33) 



Given a Galactic mass model along the line of sight, p(x), the 
prior probability of a given relative parallax is proportional to 



P(;t, vI )x J dD s Djp(D s ) J &D,D l p(D l )5[^ v 



AU AU 
(34) 

wh ere Di and D s are the l ens an d source distances. We adopt 
the lHan &Gouldl (119951 120031) model, and in Figure g] we 
plot 7r re i(M)P[7r re i(M)] versus logM to display the constraint 
placed on the mass by the measurement of 9e- The full-width 
half-maximum range is, 

log(M/M Q )=-0.7±0.4 (FWHM). (35) 

In the absence of such a measurement, the only constraint 
comes from the measurement of £e, and this is extremely 
weak, having a full width half maximum of a factor ~ 100. 
See figure 1 from Gould (2000). Indeed, as shown in that fig- 
ure, the mere supposition that the lens is a star places stronger 
constraints on the lens mass than does the measurement of t^. 

The measurements of 9% and ?e yield p K \ (eq. Ell. Since the 
distribution of fi le \ varies with D\ and D s , one could in princi- 
ple use its determination (eq. 1321 ) to place further constraints 
on combinations of these parameters and so (through eq. 1331 ) 
on the mass. In practice, for bulge sources, the distribution of 
/i le i hardly varies as a function of lens position, even when 
one considers bulge versus disk lenses. Moreover, the actual 
measured value of /j, re \ is near the peak of that distribution. 
Hence, we do not incorporate this constraint. 

The fi re \ measurement does effectively rule out a foreground 
disk source. (Without this constraint, i.e. from the CMD 
alone, the source could plausibly be a disk clump giant at 
D s ~ 5 kpc). However, for disk-disk events along this line 
of sight, the observer, lens, and source all share the same 
transverse motion due to the flat rotation curve of the Galaxy. 
Hence, only their peculiar motions relative to this rotation en- 
ter /ijei, and these are only of order 10's of kms" 1 . Hence, 
/x re i would be only a few kms" 1 kpc -1 , much slower than the 
measu red value. 

The ISumi et alJ J200 3) proper-motion measurement of the 
source independently rules out a foreground disk lens, since 
the source is moving roughly opposite to the direction of 
Galactic rotation at about p ~ —v c /Ro, where v c ~ 220 kms" 1 
and Ro = 8 kpc. In fact, this measurement by itself would be 
consistent with the source lying in the background disk, be- 
hind the bulge. However, such a scenario is virtually ruled out 



by the CMD (see Fig.|2j, which shows the source lying in or 
slightly above the bulge giant branch. If the source lay at, say, 
10 kpc, it would intrinsically be ~0.5 magnitude brighter still. 

Combinin g our m e asurem ent p le \ = 5.6 masyr" 1 , 
with the ISumi et all d2003l) measurement (p a ,^s) = 
(—2.9,-12.3) mas yr" 1 , we can effectively rule out a disk lens. 
These measurements imply |/jJ = |/fs + /-* re il ^ 7 masyr" 1 , 
whereas a disk lens would be expected to have roughly zero 
proper motion. 

5.2. Lens Luminosity M j i 

The measurement of the unlensed background flux, Fb, 
gives an upper limit to the flux from the lens. The measured 
background flux is a function of observatory and filter, and 
tends to grow with larger mean seeing. Hence, the best con- 
straint is expected to come from the observatory with the best 
seeing. In our case, this is OGLE. The OGLE Fj, is also by 
far the best constrained, in part because of the large number 
of baseline points. The OGLE background-to-source flux ra- 
tio is Ft,/F s = 0.003 ±0.009, which yields a 3 a lower limit on 
the magnitude difference of the lens and source, // — 7, > 3.6. 
For this limit to be at all relevant, the lens must be close 
to the turnoff or brighter, implying that it is close to a solar 
mass. Then, from equation (1331 . the source and lens must be 
nearly the same distance. This implies in turn that the above 
limit on apparent-magnitude difference translates directly into 
a limit on absolute-magnitude difference. Since the source 
is about 1 mag brighter than the clump, the constraint yields 
only Mi j > 2.4, which is of very limited value. 

6. MICROLENS PARALLAX tv e 

The detection of microlens parallax is marginal. We there- 
fore begin by investigating whether its tentatively detected 
value is consistent with what else is known about the lens. 
Given this orientation, and for simplicity of exposition, we 
ignore the very large error in the measurement. Only one 
component of the parallax is measured. We therefore actually 
have a limit, not a measurement, of tte > |tte,|| I = 0.86. To- 
gether with equations (j4ji and Q2> . this implies M < 0.03 M Q 
and 7r re i > 170 /ias. The source distance cannot be much more 
than D s ~ 10 kpc, partly because of the low density of stars at 
greater distances and partly because it would lie in an unpop- 
ulated portion of the CMD. Hence, 717 = 7r re i + 7r, > 270 /ias, 
i.e. D\ < 3.7 kpc. That is, the lens would be a disk brown 
dwarf. 

Apart from the small peculiar velocity of each, the lens and 
Earth are both rotating about the Galactic center at the same 
speed. Hence, the lens should be seen moving against the 
bulge at about ~ 220kms _1 towards Cygnus, which is to say 
at a position angle roughly 30° east of north. Because the 
dispersion of bulge stars is about 90 kms" 1 , this should also 
be approximately the direction of lens motion relative to the 
source. 

Since only one component of 7Te is measured, all we can 
test is the sign of this prediction. From the post-peak residu- 
als to the fit without parallax (Fig.[3J, the Earth is accelerating 
in the direction of the lens motion (thus slowing down the end 
of the event). On July 19 (roughly one month after oppo- 
sition), this is basically opposite the direction of the Earth's 
motion, and so is basically toward the west. Since the field is 
south of the ecliptic, there is also a small component of this 
(projected) acceleration toward the south. Hence, the posi- 
tion angle of the projected acceleration vector is about 260°, 
which is misaligned with the expected direction of the lens 
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motion by about 130°. That is, the expected sign of ir E u is 
opposite to what is expected. 

While it remains possible that the peculiar velocities of the 
lens and source conspire to produce this result, the statistical 
significance of the parallax measurement is not high enough 
to warrant its acceptance in the face of this strong contrary 
expectation. 

Moreover, there are at least two other possible explana- 
tions for this asymmetry apart from statistical fluctuations. 
The first is xallarap, distortions in the light curve due to ac- 
celerated motion of the source rather than the lens. Indeed, 
Smith . Mao & P aczynski (2003) showed that any parallax ef- 
fect could be mimicked by the orbital motion of the source 
around an unseen companion. When both components of 7Te 
are well measured, this possibility can be largely discounted 
because the probability of the source being in a binary with 
the same inclination, phase, and period as the Earth's orbit is 
extremely small. However, in the present case, in which all 
that is detected is a single component of acceleration, there 
is a very wide class of source binaries that could mimic the 
obs erved parallax signal. Moreover, by the arguments given 
in § 15.21 any source companion on the main-sequence would 
be undetectable in the lightcurve (other than through its ef- 
fect accelerating the source). The xallarap hypothesis could 
be checked by radial-velocity measurements. 



Still another possib le source of the asymmetry is a very 
weak binary lens. Gau di et all (|2002) detected a similarly 
weak asymmetry in OGLE-1999-BLG-36 and were able to 
model this either with parallax or with a low-mass companion 
to the lens. Thus, asymmetric residuals can be attributed to 
several effects including parallax, xallarap, and binary lenses. 
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